A Non-Overlapping Domain Decomposition Method for the Exterior Helmholtz Problem

نویسندگان

  • Armel de La Bourdonnaye
  • Charbel Farhat
  • Antonini Macedo
  • Frédéric Magoulès
  • François-Xavier Roux
چکیده

In this paper, we first show that the domain decomposition methods that are usually efficient for solving elliptic problems typically fail when applied to acoustics problems. Next, we present an alternative domain decomposition algorithm that is better suited for the exterior Helmholtz problem. We describe it in a formalism that can use either one or two Lagrange multiplier fields for solving the corresponding interface problem by a Krylov method. In order to improve convergence and ensure scalability with respect the number of subdomains, we propose two complementary preconditioning techniques. The first preconditioner is based on a spectral analysis of the resulting interface operator and targets the high frequency components of the error. The second preconditioner is based on a coarsening technique, employs plane waves, and addresses the low frequency components of the error. Finally, we show numerically that, using both preconditioners, the convergence rate of the proposed domain decomposition method is quasi independent of the number of elements in the mesh, the number of subdomains, and depends only weakly on the wavenumber, which makes this method uniquely suitable for solving large-scale high frequency exterior acoustics problems. Acoustic wave propagation problems lead to linear systems that become very large in the high frequency regime. Indeed, for most discretization methods, the mesh size h is typically chosen as one tenth of the wavelength in order to ensure a basic approximation of the physical phenomena. For this reason, many iterative solvers have been and continue to be developed for the Helmholtz problem. In this paper, we consider a domain decomposition based iterative algorithm, because of the success encountered by such methods for the solution of elliptic problems, and because they can be easily implemented on parallel computers.

برای دانلود رایگان متن کامل این مقاله و بیش از 32 میلیون مقاله دیگر ابتدا ثبت نام کنید

ثبت نام

اگر عضو سایت هستید لطفا وارد حساب کاربری خود شوید

منابع مشابه

Domain Decomposition Algorithms for an Indefinite Hypersingular Integral Equation in Three Dimensions

In this paper we report on a non-overlapping and an overlapping domain decomposition method as preconditioners for the boundary element approximation of an indefinite hypersingular integral equation on a surface. The equation arises from an integral reformulation of the Neumann screen problem with the Helmholtz equation in the exterior of a screen in R. It is well-known that the linear algebrai...

متن کامل

Restricted overlapping balancing domain decomposition methods and restricted coarse problems for the Helmholtz problem

Overlapping balancing domain decomposition methods and their combination with restricted additive Schwarz methods are proposed for the Helmholtz equation. These new methods also extend previous work on non-overlapping balancing domain decomposition methods toward simplifying their coarse problems and local solvers. They also extend restricted Schwarz methods, originally designed to overlapping ...

متن کامل

A quasi-optimal non-overlapping domain decomposition algorithm for the Helmholtz equation

This paper presents a new non-overlapping domain decomposition method for the Helmholtz equation, whose effective convergence is quasi-optimal. These improved properties result from a combination of an appropriate choice of transmission conditions and a suitable approximation of the Dirichlet to Neumann operator. A convergence theorem of the algorithm is established and numerical results valida...

متن کامل

Quasi-Optimal Convergence of Non-Overlapping Domain Decomposition Algorithm for the Helmholtz Equation

This paper presents a new non-overlapping domain decomposition method for the Helmholtz equation, whose effective convergence is quasi-optimal. These improved properties result from a combination of an appropriate choice of transmission conditions and a suitable approximation of the Dirichlet to Neumann operator. A convergence theorem of the algorithm is established and numerical results valida...

متن کامل

ذخیره در منابع من


  با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید

عنوان ژورنال:

دوره   شماره 

صفحات  -

تاریخ انتشار 1998